Broadly neutralizing antibodies against emerging delta-coronaviruses (preprint)

Porcine deltacoronavirus (PDCoV) spillovers were recently detected in children with acute undifferentiated febrile illness, underscoring recurrent zoonoses of divergent coronaviruses. To date, no vaccines or specific therapeutics are approved for use in humans against PDCoV. To prepare for possible future PDCoV epidemics, we isolated human spike (S)-directed monoclonal antibodies from transgenic mice and found that two of them, designated PD33 and PD41, broadly neutralized a panel of PDCoV variants. Cryo-electron microscopy structures of PD33 and PD41 in complex with the PDCoV receptor-binding domain and S ectodomain trimer provide a blueprint of the epitopes recognized by these mAbs, rationalizing their broad inhibitory activity. We show that both mAbs inhibit PDCoV by competitively interfering with host APN binding to the PDCoV receptor-binding loops, explaining the mechanism of viral neutralization. PD33 and PD41 are candidates for clinical advancement, which could be stockpiled to prepare for possible future PDCoV outbreaks.

Broad receptor tropism and immunogenicity of a clade 3 sarbecovirus (preprint)

Although Rhinolophus bats harbor diverse clade 3 sarbecoviruses, the structural determinants of receptor tropism along with the antigenicity of their spike (S) glycoproteins remain uncharacterized. Here, we show that the African Rinolophus bat clade 3 sarbecovirus PRD-0038 S has a broad ACE2 usage and that RBD mutations further expand receptor promiscuity and enable human ACE2 utilization. We determined a cryoEM structure of the PRD-0038 RBD bound to R. alcyone ACE2, explaining receptor tropism and highlighting differences with SARS-CoV-1 and SARS-CoV-2. Characterization of PRD-0038 S using cryoEM and monoclonal antibody reactivity revealed its distinct antigenicity relative to SARS-CoV-2 and identified PRD-0038 cross-neutralizing antibodies for pandemic preparedness. PRD-0038 S vaccination elicited greater titers of antibodies cross-reacting with vaccine-mismatched clade 2 and clade 1a sarbecoviruses compared to SARS-CoV-2 S due to broader antigenic targeting, motivating the inclusion of clade 3 antigens in next-generation vaccines for enhanced resilience to viral evolution.

Deep mutational scans of XBB.1.5 and BQ.1.1 reveal ongoing epistatic drift during SARS-CoV-2 evolution (preprint)

Substitutions that fix between SARS-CoV-2 variants can transform the mutational landscape of future evolution via epistasis. For example, large epistatic shifts in mutational effects caused by N501Y underlied the original emergence of Omicron variants, but whether such large epistatic saltations continue to define ongoing SARS-CoV-2 evolution remains unclear. We conducted deep mutational scans to measure the impacts of all single amino acid mutations and single-codon deletions in the spike receptor-binding domain (RBD) on ACE2-binding affinity and protein expression in the recent Omicron BQ.1.1 and XBB.1.5 variants, and we compared mutational patterns to earlier viral strains that we have previously profiled. As with previous RBD deep mutational scans, we find many mutations that are tolerated or even enhance binding to ACE2 receptor. The tolerance of sites to single-codon deletion largely conforms with tolerance to amino acid mutation. Though deletions in the RBD have not yet been seen in dominant lineages, we observe many tolerated deletions including at positions that exhibit indel variation across broader sarbecovirus evolution and in emerging SARS-CoV-2 variants of interest, most notably the well-tolerated {Delta}483 deletion in BA.2.86. The substitutions that distinguish recent viral variants have not induced as dramatic of epistatic perturbations as N501Y, but we identify ongoing epistatic drift in SARS-CoV-2 variants, including interaction between R493Q reversions and mutations at positions 453, 455, and 456, including mutations like F456L that define the newly emerging EG.5 lineage. Our results highlight ongoing drift in the effects of mutations due to epistasis, which may continue to direct SARS-CoV-2 evolution into new regions of sequence space.

Identification of broad, potent antibodies to functionally constrained regions of SARS-CoV-2 spike following a breakthrough infection (preprint)

The antiviral benefit of antibodies can be compromised by viral escape especially for rapidly evolving viruses. Therefore, durable, effective antibodies must be both broad and potent to counter newly emerging, diverse strains. Discovery of such antibodies is critically important for SARS-CoV-2 as the global emergence of new variants of concern (VOC) has compromised the efficacy of therapeutic antibodies and vaccines. We describe a collection of broad and potent neutralizing monoclonal antibodies (mAbs) isolated from an individual who experienced a breakthrough infection with the Delta VOC. Four mAbs potently neutralize the Wuhan-Hu-1 vaccine strain, the Delta VOC, and also retain potency against the Omicron VOCs, including recently circulating BA.4/BA.5, in both pseudovirus-based and live virus assays, and one also potently neutralizes SARS-CoV-1. The potency of these mAbs was greater against Omicron VOCs than all but one of the mAbs that had been approved for therapeutic applications. The mAbs target distinct epitopes on the spike glycoprotein, three in the receptor binding domain (RBD) and one in an invariant region downstream of the RBD in subdomain 1 (SD1). The escape pathways we defined at single amino acid resolution with deep mutational scanning show they target conserved, functionally constrained regions of the glycoprotein, suggesting escape could incur a fitness cost. Overall, these mAbs are novel in their breadth across VOCs, their epitope specificity, and include a highly potent mAb targeting a rare epitope outside of the RBD in SD1.

Evolution of antibody immunity following Omicron BA.1 breakthrough infection (preprint)

Understanding the evolution of antibody immunity following heterologous SAR-CoV-2 breakthrough infection will inform the development of next-generation vaccines. Here, we tracked SARS-CoV-2 receptor binding domain (RBD)-specific antibody responses up to six months following Omicron BA.1 breakthrough infection in mRNA-vaccinated individuals. Cross-reactive serum neutralizing antibody and memory B cell (MBC) responses declined by two- to four-fold through the study period. Breakthrough infection elicited minimal de novo Omicron-specific B cell responses but drove affinity maturation of pre-existing cross-reactive MBCs toward BA.1. Public clones dominated the neutralizing antibody response at both early and late time points, and their escape mutation profiles predicted newly emergent Omicron sublineages. The results demonstrate that heterologous SARS-CoV-2 variant exposure drives the evolution of B cell memory and suggest that convergent neutralizing antibody responses continue to shape viral evolution.

Deep mutational scans for ACE2 binding, RBD expression, and antibody escape in the SARS-CoV-2 Omicron BA.1 and BA.2 receptor-binding domains (preprint)

SARS-CoV-2 continues to acquire mutations in the spike receptor-binding domain (RBD) that impact ACE2 receptor binding, folding stability, and antibody recognition. Deep mutational scanning prospectively characterizes the impacts of mutations on these biochemical properties, enabling rapid assessment of new mutations seen during viral surveillance. However, the effects of mutations can change as the virus evolves, requiring updated deep mutational scans. We determined the impacts of all amino acid mutations in the Omicron BA.1 and BA.2 RBDs on ACE2-binding affinity, RBD folding, and escape from binding by the LY-CoV1404 (bebtelovimab) monoclonal antibody. The effects of some mutations in Omicron RBDs differ from those measured in the ancestral Wuhan-Hu-1 background. These epistatic shifts largely resemble those previously seen in the Beta variant due to the convergent epistatically modifying N501Y substitution. However, Omicron variants show additional lineage-specific shifts, including examples of the epistatic phenomenon of entrenchment that causes the Q498R and N501Y substitutions present in Omicron to be more favorable in that background than in earlier viral strains. In contrast, the Omicron substitution Q493R exhibits no sign of entrenchment, with the derived state, R493, being as unfavorable for ACE2 binding in Omicron RBDs as in Wuhan-Hu- 1. Likely for this reason, the R493Q reversion has occurred in Omicron sub-variants including BA.4/BA.5 and BA.2.75, where the affinity buffer from R493Q reversion may potentiate concurrent antigenic change. Consistent with prior studies, we find that Omicron RBDs have reduced expression, and identify candidate stabilizing mutations that ameliorate this deficit. Last, our maps highlight a broadening of the sites of escape from LY-CoV1404 antibody binding in BA.1 and BA.2 compared to the ancestral Wuhan-Hu-1 background. These BA.1 and BA.2 deep mutational scanning datasets identify shifts in the RBD mutational landscape and inform ongoing efforts in viral surveillance.

The landscape of antibody binding affinity in SARS-CoV-2 Omicron BA.1 evolution (preprint)

The Omicron BA.1 variant of SARS-CoV-2 escapes convalescent sera and monoclonal antibodies that are effective against earlier strains of the virus. This immune evasion is largely a consequence of mutations in the BA.1 receptor binding domain (RBD), the major antigenic target of SARS-CoV-2. Previous studies have identified several key RBD mutations leading to escape from most antibodies. However, little is known about how these escape mutations interact with each other and with other mutations in the RBD. Here, we systematically map these interactions by measuring the binding affinity of all possible combinations of these 15 RBD mutations (215 = 32,768 genotypes) to four monoclonal antibodies (LY-CoV016, LY-CoV555, REGN10987, and S309) with distinct epitopes. We find that BA.1 can lose affinity to diverse antibodies by acquiring a few large-effect mutations and can reduce affinity to others through several small-effect mutations. However, our results also reveal alternative pathways to antibody escape that do not include every large-effect mutation. Moreover, epistatic interactions are shown to constrain affinity decline in S309 but only modestly shape the affinity landscapes of other antibodies. Together with previous work on the ACE2 affinity landscape, our results suggest that escape of each antibody is mediated by distinct groups of mutations, whose deleterious effects on ACE2 affinity are compensated by another distinct group of mutations (most notably Q498R and N501Y).

Molecular fate-mapping of serum antibodies reveals the effects of antigenic imprinting on repeated immunization (preprint)

The ability of serum antibody to protect against pathogens arises from the interplay of antigen-specific B cell clones of different affinities and fine specificities. These cellular dynamics are ultimately responsible for serum-level phenomena such as antibody imprinting or "Original Antigenic Sin" (OAS), a proposed propensity of the immune system to rely repeatedly on the first cohort of B cells that responded to a stimulus upon exposure to related antigens. Imprinting/OAS is thought to pose a barrier to vaccination against rapidly evolving viruses such as influenza and SARS-CoV-2. Precise measurement of the extent to which imprinting/OAS inhibits the recruitment of new B cell clones by boosting is challenging because cellular and temporal origins cannot readily be assigned to antibodies in circulation. Thus, the extent to which imprinting/OAS impacts the induction of new responses in various settings remains unclear. To address this, we developed a "molecular fate-mapping" approach in which serum antibodies derived from specific cohorts of B cells can be differentially detected. We show that, upon sequential homologous boosting, the serum antibody response strongly favors reuse of the first cohort of B cell clones over the recruitment of new, naive-derived B cells. This "primary addiction" decreases as a function of antigenic distance, allowing secondary immunization with divergent influenza virus or SARS-CoV-2 glycoproteins to overcome imprinting/OAS by targeting novel epitopes absent from the priming variant. Our findings have implications for the understanding of imprinting/OAS, and for the design and testing of vaccines aimed at eliciting antibodies to evolving antigens.

Compensatory epistasis maintains ACE2 affinity in SARS-CoV-2 Omicron BA.1 (preprint)

The Omicron BA.1 variant emerged in late 2021 and quickly spread across the world. Compared to the ancestral Wuhan Hu-1 strain and other pre-Omicron SARS-CoV-2 variants, BA.1 has many mutations, a number of which are known to enable antibody escape. Many of these antibody-escape mutations individually decrease the spike receptor-binding domain (RBD) affinity for ACE2 in the background of early SARS-CoV-2 variants, but BA.1 still binds ACE2 with high affinity. The fitness and evolution of the BA.1 lineage is therefore driven by the combined effects of numerous mutations. Here, we systematically map the epistatic interactions between the 15 mutations in the RBD of BA.1 relative to the Wuhan Hu-1 strain. Specifically, we measure the ACE2 affinity of all possible combinations of these 15 mutations (215 = 32,768 genotypes), spanning all possible evolutionary intermediates from the ancestral Wuhan Hu-1 strain to BA.1. We find that immune escape mutations in BA.1 individually reduce ACE2 affinity but are compensated by epistatic interactions with other affinity-enhancing mutations, including Q498R and N501Y. Thus, the ability of BA.1 to evade immunity while maintaining ACE2 affinity is contingent on acquiring multiple interacting mutations. Our results implicate compensatory epistasis as a key factor driving substantial evolutionary change for SARS-CoV-2 and are consistent with Omicron BA.1 arising from a chronic infection.

Mosaic RBD nanoparticles protect against multiple sarbecovirus challenges in animal models (preprint)

To combat future SARS-CoV-2 variants and spillovers of SARS-like betacoronaviruses (sarbecoviruses) threatening global health, we designed mosaic nanoparticles presenting randomly-arranged sarbecovirus spike receptor-binding domains (RBDs) to elicit antibodies against conserved/relatively-occluded, rather than variable/immunodominant/exposed, epitopes. We compared immune responses elicited by mosaic-8 (SARS-CoV-2 and seven animal sarbecoviruses) and homotypic (only SARS-CoV-2) RBD-nanoparticles in mice and macaques, observing stronger responses elicited by mosaic-8 to mismatched (not on nanoparticles) strains including SARS-CoV and animal sarbecoviruses. Mosaic-8 immunization showed equivalent neutralization of SARS-CoV-2 variants including Omicron and protected from SARS-CoV-2 and SARS-CoV challenges, whereas homotypic SARS-CoV-2 immunization protected only from SARS-CoV-2 challenge. Epitope mapping demonstrated increased targeting of conserved epitopes after mosaic-8 immunization. Together, these results suggest mosaic-8 RBD-nanoparticles could protect against SARS-CoV-2 variants and future sarbecovirus spillovers.

The SARS-CoV-2 Delta variant induces an antibody response largely focused on class 1 and 2 antibody epitopes (preprint)

Exposure histories to SARS-CoV-2 variants and vaccinations will shape the specificity of antibody responses. To understand the specificity of Delta-elicited antibody immunity, we characterize the polyclonal antibody response elicited by primary or mRNA vaccine-breakthrough Delta infections. Both types of infection elicit a neutralizing antibody response focused heavily on the receptor-binding domain (RBD). We use deep mutational scanning to show that mutations to the RBD's class 1 and class 2 epitopes, including sites 417, 478, and 484-486 often reduce binding of these Delta-elicited antibodies. The anti-Delta antibody response is more similar to that elicited by early 2020 viruses than the Beta variant, with mutations to the class 1 and 2, but not class 3 epitopes, having the largest effects on polyclonal antibody binding. In addition, mutations to the class 1 epitope (e.g., K417N) tend to have larger effects on antibody binding and neutralization in the Delta spike than in the D614G spike, both for vaccine- and Delta-infection-elicited antibodies. These results help elucidate how the antigenic impacts of SARS-CoV-2 mutations depend on exposure history.

Shifting mutational constraints in the SARS-CoV-2 receptor-binding domain during viral evolution (preprint)

SARS-CoV-2 has evolved variants with substitutions in the spike receptor-binding domain (RBD) that impact its affinity for ACE2 receptor and recognition by antibodies. These substitutions could also shape future evolution by modulating the effects of mutations at other sites--a phenomenon called epistasis. To investigate this possibility, we performed deep mutational scans to measure the effects on ACE2 binding of all single amino-acid mutations in the Wuhan-Hu-1, Alpha, Beta, Delta, and Eta variant RBDs. Some substitutions, most prominently N501Y, cause epistatic shifts in the effects of mutations at other sites, thereby shaping subsequent evolutionary change. These epistatic shifts occur despite high conservation of the overall RBD structure. Our data shed light on RBD sequence-function relationships and facilitate interpretation of ongoing SARS-CoV-2 evolution.

Selection analysis identifies unusual clustered mutational changes in Omicron lineage BA.1 that likely impact Spike function (preprint)

Among the 30 non-synonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (i) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (ii) interactions of Spike with ACE2 receptors, and (iii) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron over all previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

An antibody-escape calculator for mutations to the SARS-CoV-2 receptor-binding domain (preprint)

A key goal of SARS-CoV-2 surveillance is to rapidly identify viral variants with mutations that reduce neutralization by polyclonal antibodies elicited by vaccination or infection. Unfortunately, direct experimental characterization of new viral variants lags their sequence-based identification. Here we help address this challenge by aggregating deep mutational scanning data into an "escape calculator" that estimates the antigenic effects of arbitrary combinations of mutations to the viruss spike receptor-binding domain (RBD). The calculator can be used to intuitively visualize how mutations impact polyclonal antibody recognition, and score the expected antigenic effect of combinations of mutations. These scores correlate with neutralization assays performed on SARS-CoV-2 variants, and emphasize the ominous antigenic properties of the recently described Omicron variant. An interactive version of the calculator is at https://jbloomlab.github.io/SARS2_RBD_Ab_escape_maps/escape-calc/, and we provide a Python module for batch processing.

Antibody-mediated broad sarbecovirus neutralization through ACE2 molecular mimicry (preprint)

Understanding broadly neutralizing sarbecovirus antibody responses is key to developing countermeasures effective against SARS-CoV-2 variants and future spillovers of other sarbecoviruses. Here we describe the isolation and characterization of a human monoclonal antibody, designated S2K146, broadly neutralizing viruses belonging to all three sarbecovirus clades known to utilize ACE2 as entry receptor and protecting therapeutically against SARS-CoV-2 beta challenge in hamsters. Structural and functional studies show that most of the S2K146 epitope residues are shared with the ACE2 binding site and that the antibody inhibits receptor attachment competitively. Viral passaging experiments underscore an unusually high barrier for emergence of escape mutants making it an ideal candidate for clinical development. These findings unveil a key site of vulnerability for the development of a next generation of vaccines eliciting broad sarbecovirus immunity.

A SARS-CoV-2 variant elicits an antibody response with a shifted immunodominance hierarchy (preprint)

Many SARS-CoV-2 variants have mutations at key sites targeted by antibodies. However, it is unknown if antibodies elicited by infection with these variants target the same or different regions of the viral spike as antibodies elicited by earlier viral isolates. Here we compare the specificities of polyclonal antibodies produced by humans infected with early 2020 isolates versus the B.1.351 variant of concern (also known as Beta or 20H/501Y.V2), which contains mutations in multiple key spike epitopes. The serum neutralizing activity of antibodies elicited by infection with both early 2020 viruses and B.1.351 is heavily focused on the spike receptor-binding domain (RBD). However, within the RBD, B.1.351-elicited antibodies are more focused on the "class 3" epitope spanning sites 443 to 452, and neutralization by these antibodies is notably less affected by mutations at residue 484. Our results show that SARS-CoV-2 variants can elicit polyclonal antibodies with different immunodominance hierarchies.

ACE2 binding is an ancestral and evolvable trait of sarbecoviruses (preprint)

Two different sarbecoviruses have caused major human outbreaks in the last two decades. Both these sarbecoviruses, SARS-CoV-1 and SARS-CoV-2, engage ACE2 via the spike receptor-binding domain (RBD). However, binding to ACE2 orthologs from humans, bats, and other species has been observed only sporadically among the broader diversity of bat sarbecoviruses. Here, we use high-throughput assays to trace the evolutionary history of ACE2 binding across a diverse range of sarbecoviruses and ACE2 orthologs. We find that ACE2 binding is an ancestral trait of sarbecovirus RBDs that has subsequently been lost in some clades. Furthermore, we demonstrate for the first time that bat sarbecoviruses from outside Asia can bind ACE2. In addition, ACE2 binding is highly evolvable: for many sarbecovirus RBDs there are single amino-acid mutations that enable binding to new ACE2 orthologs. However, the effects of individual mutations can differ markedly between viruses, as illustrated by the N501Y mutation which enhances human ACE2 binding affinity within several SARS-CoV-2 variants of concern but severely dampens it for SARS-CoV-1. Our results point to the deep ancestral origin and evolutionary plasticity of ACE2 binding, broadening consideration of the range of sarbecoviruses with spillover potential.

Multivalent designed proteins protect against SARS-CoV-2 variants of concern (preprint)

Escape variants of SARS-CoV-2 are threatening to prolong the COVID-19 pandemic. To address this challenge, we developed multivalent protein-based minibinders as potential prophylactic and therapeutic agents. Homotrimers of single minibinders and fusions of three distinct minibinders were designed to geometrically match the SARS-CoV-2 spike (S) trimer architecture and were optimized by cell-free expression and found to exhibit virtually no measurable dissociation upon binding. Cryo-electron microscopy (cryoEM) showed that these trivalent minibinders engage all three receptor binding domains on a single S trimer. The top candidates neutralize SARS-CoV-2 variants of concern with IC50 values in the low pM range, resist viral escape, and provide protection in highly vulnerable human ACE2-expressing transgenic mice, both prophylactically and therapeutically. Our integrated workflow promises to accelerate the design of mutationally resilient therapeutics for pandemic preparedness.

Stabilization of the SARS-CoV-2 Spike receptor-binding domain using deep mutational scanning and structure-based design (preprint)

The unprecedented global demand for SARS-CoV-2 vaccines has demonstrated the need for highly effective vaccine candidates that are thermostable and amenable to large-scale manufacturing. Nanoparticle immunogens presenting the receptor-binding domain (RBD) of the SARS-CoV-2 Spike protein (S) in repetitive arrays are being advanced as second-generation vaccine candidates, as they feature robust manufacturing characteristics and have shown promising immunogenicity in preclinical models. Here, we used previously reported deep mutational scanning (DMS) data to guide the design of stabilized variants of the RBD. The selected mutations fill a cavity in the RBD that has been identified as a linoleic acid binding pocket. Screening of several designs led to the selection of two lead candidates that expressed at higher yields than the wild-type RBD. These stabilized RBDs possess enhanced thermal stability and resistance to aggregation, particularly when incorporated into an icosahedral nanoparticle immunogen that maintained its integrity and antigenicity for 28 days at 35-40{degrees}C, while corresponding immunogens displaying the wild-type RBD experienced aggregation and loss of antigenicity. The stabilized immunogens preserved the potent immunogenicity of the original nanoparticle immunogen, which is currently being evaluated in a Phase I/II clinical trial. Our findings may improve the scalability and stability of RBD-based coronavirus vaccines in any format and more generally highlight the utility of comprehensive DMS data in guiding vaccine design.

The SARS-CoV-2 mRNA-1273 vaccine elicits more RBD-focused neutralization, but with broader antibody binding within the RBD (preprint)

The emergence of SARS-CoV-2 variants with mutations in key antibody epitopes has raised concerns that antigenic evolution will erode immunity. The susceptibility of immunity to viral evolution is shaped in part by the breadth of epitopes targeted. Here we compare the specificity of antibodies elicited by the mRNA-1273 vaccine versus natural infection. The neutralizing activity of vaccine-elicited antibodies is even more focused on the spike receptor-binding domain (RBD) than for infection-elicited antibodies. However, within the RBD, binding of vaccine-elicited antibodies is more broadly distributed across epitopes than for infection-elicited antibodies. This greater binding breadth means single RBD mutations have less impact on neutralization by vaccine sera than convalescent sera. Therefore, antibody immunity acquired by different means may have differing susceptibility to erosion by viral evolution.